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We introduce a new path integral representation for slave bosons in the radial gauge which is 
valid beyond the conventional fluctuation corrections to a mean-field solution. For electronic lattice 
models, defined on the constrained Fock space with no double occupancy, all phase fiuctuations of 
the slave particles can be gauged away if the Lagrange multipliers which enforce the constraint on 
each lattice site are promoted to time-dependent fields. Consequently, only the amplitude (radial 
part) of the slave boson fields survives. It has the special property that it is equal to its square in the 
physical subspace. This renders the functional integral for the radial field Gaussian, even when non- 
local Coulomb-type interactions are included. We propose i) a continuum integral representation 
for the set-up of further approximation schemes, and ii) a discrete representation with an Ising-like 
radial variable, valid for long-ranged interactions as well. The latter scheme can be taken as a 
starting point for numerical evaluations. 
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I. INTRODUCTION 



The significance of strong electronic correlations in metals has long been acknowledged when intractable enigmas of 
solid state physics such as itinerant magnetism, the correlation-induced metal-insulator transition, unconventional 
superconductivity, or spin liquid states and non-Fermi liquid behavior of electronic excitations were embraced. The 
scope of the "strongly correlated physics" is nowadays enormous; it extends to the fractional quantum Hall effect and 
high-Tc superconductivity, and also in mesoscopic systems an increasing number of investigations take this route since 
certain scenarios with strong local interactions can be modeled more directly in experiments on a mesoscopic scale. 
The efforts to fathom the underlying mechanisms resulted only in partial success. Especially, models with strong local 
interactions of itinerant electrons resisted a thorough understanding beyond certain limiting, and often unphysical, 
cases. 

The celebrated Hubbard model epitomizes the two opposing characters of narrow band electrons: electron hopping 
(with energy scale t) supports the itinerant character, suitably expressed in momentum space, and favors a metallic 
character. In contrast to this, the local on-site Coulomb interaction U may drive the electrons, depending on band- 
filling, into a Mott insulating state, with strong magnetic correlations. Standard perturbation theory is not applicable 
in the transition regime (for a recent discussion see j^]), and whether it even was valid for any finite U in two space 
dimensions |^ was a longstanding problem. On the other hand, perturbation theory from the atomic limit, with no 
kinetic terms, is not feasible in a straightforward manner since the tremendous degeneracy for low energy states of 
the fully localized limit renders the usual perturbation theory inapplicable. Furthermore, a hopping expansion cannot 
build on Wick's theorem since the zeroth order Hamiltonian is not bilinear in the canonical electron fields, though 
progress has been achieved recently Nonetheless, for t ^ J7 a Schrieffer- Wolff transformation may be performed 

to the appropriate order in t/U. This may be visualized as a rotation of the states in Hilbert space such that transitions 
between different subspaces with fixed number of doubly occupied sites are suppressed to the corresponding order in 
t/U. Equivalently, one usually applies a canonical transformation to the Hamilton operator in order to generate an 
effective Hamiltonian of which the low energy part is explicit. The Hubbard model is reduced in this way to a 
t-J model (plus pair-hopping) up to linear order in t/U — or to a "i-model" in zeroth order. The t-model is just 
the kinetic term, however it is now a correlated hopping in the sense that double occupancies are excluded. This 
constrained hopping can only be envisioned in real space. The t-J model includes, beyond the constrained hopping, 
a Heisenberg spin-exchange and a short range density-density correlation term of scale J ~ t'^ /U. Although all these 
strong coupling models present the local spin and charge excitations in an appealing way, it is still a long way to 
extract the true (delocalized) low-energy excitations. 

In the past decades considerable effort has been invested into the construction of techniques which work in the 
constrained Fock space. Already Hubbard introduced projection operators onto local (site-) states with no electron 
(|0)), one electron with spin a (jcr)) or double occupancy (|2)). These so-called Hubbard X-operators belong to a 
graded Lie-algebra, and a canonical many-particle Green's function formalism is not applicable (for a detailed review of 
functional integral representations and a suitable choice of coherent states for X-operators see ref . Q ) . The evaluation 
scheme by Hubbard shaped the appropriate picture of lower and upper "Hubbard bands" for low and high energy 
(in-) coherent excitations, respectively. Yet it suffered from inconsistencies, e.g. the failure to satisfy the Luttinger 
theorem. 

An alternative approach, which also implements a projection onto local states, is the slave boson formalism 
Electron creation and annihilation operators are represented by composite operators whereby an electron is annihilated 
either by the creation of a hole site, accompanied by the annihilation of a singly occupied site with corresponding spin 
or the annihilation of a doubly occupied site, accompanied by the creation of a singly occupied site with opposite spin. 
Several choices exist how to ascribe fermionic and bosonic degrees of freedom to these operators, for a discussion see 
. Although the operators are now canonical fermionic and bosonic operators, they live in a constrained Fock space 
since the sum of their respective particle numbers (the pseudo charge) has to be unity for each site. Mean field or 
saddle point approximations satisfy the constraint only on the average and, furthermore, lead to Bose condensation 
of the slave bosons |l0|| . Actually, the saddle point approximation becomes exact in the unphysical limit A/" ^ oo, 
where Af is the spin degeneracy. Corrections to the mean field result were calculated in the first order of the expansion 
parameter l/M — the constraint is thereby replaced by a "soft constraint" with pseudo charge equal to A/'/2. For 
Af = oo, the bosons condense, that is, they attain a finite complex expectation value, the square of which is the 
density of empty sites. 

In the Barnes representation, only the charge degree of freedom is represented by the slave boson. In contrast, both 
spin and charge degrees of freedom are carried by bosons (and fermions) in the Kotliar and Ruckenstein representation 
while in the spin and charge rotation invariant representation they are carried by bosons only. The mean- 
field approach to the Kotliar and Ruckenstein representation has been quite successful when compared to numerical 
simulations. Indeed, the paramagnetic mean-field solution is equivalent to the Gutzwiller approximation Even 
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though this is the saddle point of an action, it turns out to obey a variational principle in the limit of large spatial 
dimensions, where the Gutzwiller approximation and the Gutzwiller wave function are identica l When the 

dimension is reduced, they are not known to differ markedly. Comparison of ground state energies [|l2,Q, and charge 
structure factors jl4j with Quantum Monte Carlo data show excellent agreement. A variety of saddle points have 
been investigated, in particular in the vicinity of the Mott transition p5|-p^. The Mott gap can be calculated as well 
and comparison with experimental data in La^jYi-ajTiOs ||2^ displays a very good agreement Nevertheless 
the computation of fluctuations, in particular in the magnetic channel, remains problematic. 

Moreover such a finite expectation value of the complex slave boson field violates Elitzur's theorem according 
to which, in a theory with local gauge invariance, only locally gauge invariant operators may acquire a non-zero 
expectation value. Correspondingly, for finite Af, the phase fluctuations of the boson field suppress the condensation. 
A zero expectation value of the bosonic field, however, is problematic as a starting point since the bosonic field has 
to connect to the number of empty sites, in the end. It was early on recognized that these difficulties are avoided in 
a polar decomposition where the zero mode of the phase fluctuations may be singled out and the expectation value 
of the radial amplitude is, in general, non-vanishing pl| , p^ (a more recent, detailed discussion of l/A/'-corrections in 
the "radial gauge" is found in the review by Arrigoni et al. ||23| ). Although subtleties are associated with the measure 
and the correct translation of the time discretized functional integral into the continuum limit version, the result in 
l/A/" is not different from that of an X-operator approach |Q. 

Despite the progress in the understanding of the fluctuations around the mean field saddle point, the physically 
most interesting case with the largest possible spin fluctuations, Af — 2, remains inaccessible. Moreover, the saddle 
point solution is generically a Fermi liquid solution and, therefore, non-Fermi liquid behavior is necessarily non- 
perturbative in 1/J\f. This motivated the Karlsruhe group to device a slave-boson formalism which builds on 
canonical many-particle techniques, is locally gauge-invariant by construction and, most importantly, allows for Fermi 
liquid and non-Fermi liquid behavior for the considered (multi-channel) Anderson impurity problem. It permits to 
recover the physical quantities over the complete temperature range, from the high temperature local moment regime 
to the low temperature correlated many-body state. In spite of this major success of slave boson theory, it is still a 
formidable problem to set up a similar diagrammatic evaluation for strongly correlated itinerant electrons. 

It is quite instructive to compare the implementation of the constraint in this diagrammatic formalism with that in 
the functional integral radial gauge representation. In the first approach, local gauge invariance is guaranteed through 
Ward identities in a conserving approximation. This enforces the constraint only in so far as transitions between sectors 
of the Fock space with different pseudo charge are suppressed during the time evolution of a considered many particle 
state. The projection onto the physical sector of the Fock space is achieved in a last step, when physical expectation 
values are calculated, using the Abrikosov procedure of sending the Lagrange multiplier of the constraint to infinity. 
This may be contrasted with the functional integral calculation in the radial gauge: here the phase is fixed and the 
amplitude of the boson will in turn strongly fiuctuate. The pseudo charge conservation will not be satisfied a priori. 
This is precisely the reason why one has to project onto the physical subspace at each time step (and, of course, at 
each site) which is accomplished by promoting the Lagrange multiplier of the constraint into a field. 

In this paper we reexamine functional integrals for strongly correlated electrons with A/" = 2 in the radial gauge. 
We restrain our considerations to "C/ — > cxd" , that is the limit of strong local interaction between electrons so that 
the number of doubly occupied sites is zero. The Hubbard model is reduced to either the t— ot t ~ J model. In this 
limit, the bosonic fields in radial representation reduce to their respective (real) amplitude since all time derivatives 
of the corresponding phases can be absorbed in a (time-dependent) Lagrange multiplier field. 

It is possible to construct the functional integral either as an integration over continuous values of the boson radial 
amplitudes and the Lagrange multiplier, or, equivalently one can implement the constraint first and thereby reduce 
the integration over the slave boson amplitude to a discrete sum over the values of "hole occupancy", {0, 1}. The first 
method is more appropriate for approximation schemes which further manipulate the path integral analytically, and 
the latter method is to be taken as a starting point for numerical calculations in the strong coupling limit which do 
not build on a Hirsch decoupling |^ . 

In general, there are no ab initio advantages of one functional integral representation over another representation. 
It depends on the approximation scheme if a concrete representation is suitable. For example, when l/A/'~fluct nations 
are considered the radial gauge consistently separates the zero-mode fluctuations. Moreover, for any A/", the square 
of the slave boson amplitude is equal to the amplitude itself since the constraint projects the amplitude onto the 
values {0, 1}. In the radial gauge, this renders the integrations over the real fields Gaussian, albeit the bosons are still 
coupled to the fermionic Grassmann fields. This projection property also ensures that the measure does not produce 
the cumbersome logarithmic terms in the action and, more significantly, non-local density-density interactions turn 
out to be Gaussian terms in the functional integral. The latter is not valid in the Cartesian gauge. These observations 
should suffice as motivation to readdress the construction of the slave boson functional integral in the radial gauge. 

In Se ctio n |l^ we introduce the correct functional integral representation in the radial gauge, in the limit U — *■ oo. In 
Section III wc propose an alternative discrete representation which may serve for numerical simulations. In Section IV 
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we extend these representations to lattice models, with (possibly) long-ranged density-density interactions. Several 
calculations for simple models are carried out to illustrate this new method. 



II. FUNCTIONAL INTEGRAL REPRESENTATION IN THE RADIAL GAUGE. 



In any slave boson representation of strongly correlated electron systems there is a gauge symmetry group as 
a consequence of the redundancy of the representation of the physical electron operator. The first slave boson 
representation is due to Barnes 0|. It was introduced to tackle the single-impurity Anderson model. The local 
physical electron operator cj. is, in the limit of infinitely strong on-site interaction, rewritten as a product of auxiliary 
fields fl and b: 

4 = bfl (1) 

Here is a doublet of canonical fermionic fields, and b a canonical bosonic field. We stay in the physical sector of 
the Fock space provided the auxiliary fields fulfill the constraint relation: 

Q = ^/t/, + 5tfe=l (2) 

(7 

In a functional integral representation of the partition sum Z, the constraint is enforced via the integration over a 
Lagrange multiplier A (for a standard reference see pTf): 



(3) 



with 

£'{T)=Cf{T)+Ct{T)+£t{T)+ £„ioc (r ) (4) 

and 

/:/(T)^^/t(r)(a,-M + *A)/.(r) 

CbiT)=b^ir){dr+i\)biT) (5) 

The last two contributions to the Lagrange density are hopping of electrons (£t) and the non-local part of the electron- 
electron interaction (£nioc), respectively. The chemical potential, which fixes the average number of electrons, is /i. 
The space dependence (lattice site index i) of the fields and of A is suppressed in Eq. (y|), for simplicity. Since the 
pseudo charge Q, defined in Eq. (^, is a conserved quantity, the action is invariant under a group of local U{1) gauge 
transformations p2| . In the contimmm limit, it reads: 

/.(r)^/.(r)-^(r)e^'^M 
5(r) ^ 6(r) = 6(T)e*'^(^' 

A ^ A(r) = A - ^(r) . (6) 

The physical electron operators are invariant under this transformation. It is well known that A has to be continued 
into the complex plane, and the integration contour has to be shifted into the lower half-plane in order to ensure 



convergence of the functional integral |28 



A ^ A - iXa (7) 

Moreover, the phase of the bosonic field may be gauged away by promoting the constraint into a field. The remaining 
degree of freedom of the boson is consequently its amplitude. Having gauged away the phase of the slave boson brings 
us to the radial gauge. 

Certainly, the expectation value of any physical observable has to be gauge invariant. In the conventional Cartesian 
representation — with the integration over both, the real and imaginary components of the bosonic field — the particle 
number of bosons (holes) is fixed within a gauge-invariant, conserving approximation p^ ]. In the radial representation, 
however, the gauge is usually fixed in order to reduce the bosonic field to a real field, and the number of bosons is 
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fluctuating. Only in the final integration over the Lagrange multiplier field A(t), the particle number conservation is 
ensured at each time step. This situation is somewhat reminiscent of the theoretical modeling of superconductivity 
where, in BCS theory, the phase of the condensate is fixed, but gauge invariance is restored in RPA Nonetheless 
the situation for superconductivity is distinct because there, a global gauge is broken, whereas here a broken local 
gauge has to be restored. 

In a specific gauge, the correctness of the representation can be verified through the evaluation of the partition 
sum, and any correlation function, in the atomic limit. In this case the calculation is straightforward. For example, 
the partition sum reads: 



J-7:/p 27r 



(8) 



where Zf = [I + e-'3(»-^+^o-A')]2 r^^^ Zf, = [1 - e"'^(*^+^o)]-i^ ]\[Qtc that the projection onto the physical subspace 
{J dX) follows from a contour integral in the complex z-plane, with z = e^^^^, and it only picks up the residue at 
z = 0. 

One may try to carry out this program in the polar representation, namely by introducing the amplitude r„ and 
the phase of the bosonic field: 



db'„db'' 



d{rl 



d(pn 

2tt 



(9) 



The bosonic part of the partition sum can be calculated ^^1; the result is again correctly Zi,. Nevertheless, in 
this way one does not promote A into a field, and the phases ip„ are notgauged away and, more importantly, the 
action to lowest order in the time step, is not a bilinear form in r and (p [ pO[ . 

We now start to set up the path integral representation in the "radial gauge". In the first place we observe that, 
contrary to speculations that the radial representation simply results from a straightforward coordinate transformation, 
it needs to be set up on a discretized time mesh from the beginning. A translation which respects the requirement 
that the integrals converge, irrespective of the sequence of the various integrations, is the following: 



lim lim lini 

N^oo e^0+ n^cx) 



D 



S ■ dXr, 
27r 



(10) 



where 



N 



n— 1 a 

N 



(11) 



n=l 



and 5 = (3/N with N the number of time steps. 

We proceed with the evaluation of the atomic partition sum, step by step, in order to convince the reader of the 
correctness of this representation: integration over the Grassmann variables [fa,f^] yields (we implicitly understand 
that the limit iV ^ cxd has to be taken at the end of the calculation): 



/ N n 
Zat = lim lim TT / 



5dX„ 



d{rl)e 



N 



1+ n (1 - '5(iA™ + Ao - /i)) 



^A-(iA„ + Ao)(<-e-l) 



{n>Xo> 0) 



(12) 



The path integral representation, and in particular the Trotter-Suzuki decomposition, only makes sense when |fa| <g; 1, 
where x denotes anything that multiplies 5 in the course of the calculation. Clearly, a finite |A„| < £7 is required in 
order to reexponentiate the fermionic contribution, or, in other words, to introduce a small parameter, allowing to 
neglect terms of order 0{5'^), for our purpose. The square bracket in Zat therefore may be written as 



Zf 



N 
n=l 



-i5(iA„+Ao— /j) 



(13) 
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and we now obtain 



Zat = lim lim 



' N 

n 



lim lim 



n 



n 



vn=l • 
„5(iA„+Ao)(l+e) 

(5(iA„ + Ao) 



SdXn 

2ti 



N g5(iA„+Ao)e 



\ 5{iXn + Ao) 



+ e 



AT 

n 



-5(i:A„ + Ao)(l-e) 
5(iA„ + Ao) 



(14) 

A motivation, 



The origin of e is now clear: it prescribes how to close the contour for the second term in Eq. mA 
how to find e and the bosonic Lagrange density in radial representation, is presented in AppendixRI. 

In each term, the existence of the integrals is guaranteed by the shift of the A-integration contour into the lower 
complex plane induced by a positive Aq. The residue of the pole at A = iAo of the contour integration in the first 
two terms is finite and yields exactly the two contributions of Zat : the hole line and the two particle lines with single 
occupancy of Fig. 1, respectively. The third term in Eq. (p^, corresponding to double occupancy, is annihilated 
through the A-integration, as expected from the implementation of the constraint. 

It is instructive to invert the sequence of the integrations. 



lim 



SdXn 
2n 



lim I I / d(r 
hmny dirl) 

e-^^"('--.-i)^(r2 - 1) + 2e-'^°^-"^-'^{rl ~ e)e/'' + er'^^^^'^+'^irl + l)e^'^ 



(15) 



where S is the Dirac (5-function. The last line is equivalent to: 

Zat = l + 2e'3^ . 

In this way it appears that, after the projection has been performed in the first step, r„ may only take two integer 
values, or 1 in the last line of Eq. (p^). We thus arrive at the conclusion that the integrations over the bosonic and 
the constraint fields can be replaced by the simple procedure demonstrated in Eq. ([l5|). 

To summarize, we have introduced a new path integral representation for slave bosons in the radial gauge. It 
is defined on a discretized time mesh, and has a well-defined continuum limit. Furthermore, we can take a more 



conventional approach (Eq. (10)), based on a coherent state functional integral for both fermions and bosons, or use 
an Ising-like variable representing the amplitude of empty sites (Eq. ([T5|)). The extension to non-local terms in the 
Lagrangian will be discussed in Section IV where we also present a technique how to handle the integration over the 
amplitude field r^. 



III. DISCRETE REPRESENTATION 



Besides the above integral representation with continuous fields, which may be used to set up various approximation 
schemes, one can alternatively obtain a discrete representation. Indeed the above calculation of the partition sum 
already suggests that it can be obtained without having to integrate over continuous values of r„. Moreover integration 
over A„ mostly has the effect of "picking up" the appropriate value of r„. This is achieved more elegantly by 
extending Coleman's projection scheme pl| to the case of a time-dependent constraint. For this purpose one introduces 
^„ = e~'^(*^"+^''), and rewrites the partition sum in the following compact form: 

(16) 

€i=...=?jv=o 
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with 



N 



^/ = E E -^L [fn,a - fn-lAne'''] (17) 
n—1 fT 

In Eq. (|6|) and all following relations, J2{r =o 1} denotes the "path integral" over the field r which now is reduced 
to the set of discrete values r„ = 0, 1 for n = 1, . . . ,N. We easily check that the atomic limit is reproduced: 



Nr. \ / N 

a 



{r„=0,l}„ \n=l 



= ...=4iv=0 

= 1 + 26^^^ . (18) 

This discretized representation is certainly the most direct way to construct the partition sum — and later on the 
Green's function — as a sum of free fermionic paths, controlled by a set of Ising variables ri to r^. 

A. Representation of the Grassmann fields for a physical electron. 

In the radial gauge we have the following assignment of electron fields to the previously introduced auxiliary fields: 

= r ft 

Cn,a — '^71+1 f 71, a (1^) 

The product composition is just the expression Eq. (^ — only the choice of the time steps has to be confirmed. 
This is easily derived from the requirement that in order to attain a non-zero result for the Green's function we 
need "complete chains" in the Grassmann integration. We first clarify the term "complete chain" in the following 
paragraph, and then convince ourselves in Subsection IIIB that the choice of time steps in Eq. (|l9|) is necessary to 



obtain a finite electron propagator. 

"Complete chain" denotes a product of Grassmann numbers over all time steps, either of the form J^^^j^ /rLo-/,! o- 

or, equivalently, as Y\n=i fn afn-i,cr which both result from the expansion of the exponentiated action, and only 
the integration of such complete chains results in a non-zero value of the path integral. The first form can be 
characterized as a "hole chain" since it corresponds to the hole line in Fig. 1 (top line). Actually, this hole contribution 
to Zat with value 1 is formed by the product of two hole chains (0^=1 T)(n^=i fndfl i)- presented 

slave boson scheme, this double hole chain is matched by a complete chain of r„ variables with value 1 (the first 
contribution in the last line of Eq. (p^). The second arrangement 0^=1 /n cr/n-i,o- then should be interpreted 
as a "particle chain" which carries spin a. Indeed, the dashed lines with value exp(/i/3) in Fig. 1 are formed by 
(n^=i /n,-c/n.-cr)(n^=i /n.cr/n-i.o-e'''') with (T —1,1 for the two lines, respectively. Finally, the third contribution 
in the last line of Eq. ( [T5[ ) is generated from a double particle fine (11^=1 /ri,-o-/n-i,-CTe'^^)(nn=i /n.^/n-i.^e*^) but 
it is projected out via S{r'^ + 1), as we work in the constrained Fock space with no double occupancies. 

B. Calculation of the Green's function. 

For the Green's function Gc{Tf — Ti) — with electron creation at time step and annihilation at t/ — the 
complete chains again have to exist. However they are built from a distinct arrangement of Grassmann variables from 
the Lagrangian and the physical electron creation and annihilation operators. In this case, the hole chain runs only 
from time step to the time step when the cr-electron is created (visualized in Fig. 2, where the electron creation is 
at time step 1). The u-hole chain is then replaced by a cr-particle chain in between time steps 1 and m. This product 
of hole- and particle chains still results in a finite value for the Grassmann integrations since the electron creation 
and annihilation operators at time steps 1 and m provide the missing Grassmann numbers in order to produce a 
"complete chain" from to A^. 

In order to verify this arrangement and to understand the shift of the time arguments by one time step in the 



relation for Cn,a in Eq-(19), we explicitly carry out the evaluation of the Green's function 

-ZGa{m - 1) = (c,„,o-4,o-) = {^777+1 fni.arifl^) (20) 



6 



in the atomic limit. Obviously the cr-fermion exists for m — 1 time steps. When we integrate out the Grassmann fields 
as above for Zat we have: 



-^atGat,.(m-l)= J2 (l[l^^n]rm+iri(l+l[C„,eA D-\ 

{r„=0,l}„ \n=l / \ m=l / 



(21) 



Here the integration over the —cr-fermion produces one factor (1 + 11™=! (.me-^'^) = det [D], and the integration over 
the CT-fermion yields the propagator matrix- element D^^i multiplied by the second factor det [D]. This propagator is: 

m 

D-,],^l[{^.e"^)/det[D] (22) 

1=2 

Next we calculate the derivatives at each time step and obtain 

-^atGat,.(m - 1) = J2 '^™+ir-i((5.i,i<5.,,o...'^r„,o5r„+i,i...<5r„,i)e*^^'"~'^ (23) 

{r„=04}n 

As anticipated above, we produced a chain of r-values, with = for 2 < i < m. We find the correct expression for 
this projected atomic Green's function 

-^atGat,.(m - 1) = e^("-i)* . (24) 



If we had erroneously taken Cn.a — 'Tnfn.a in Eqs- (|19| ) and (|2Cl|), we would have to replace r,„+i by r,,,, in Eq. (|23|), 
and the result would be incorrectly zero, due to the <5r„,o in the r— chain. Higher order correlation functions can be 
calculated in a similar way. 

We summarize that the expressions for the decomposition of the c— operators (Eq. (p^)), together with the expression 
of the local part of the action (Eq. (^l])) and the measure of Eq. (|o|) form a complete representation of the slave 
boson functional integral in the radial gauge in the limit of infinitely strong local interaction. In the following chapter 
we will discuss non-local terms in the Lagrangian, using exactly this representation. 



IV. NON-LOCAL LAGRANGIAN 



Two types of non-local terms have to be included for a full discussion of physical models for correlated electrons: 
kinetic terms, i.e. hopping of electrons, and non-local interactions. 
We first consider kinetic terms of the type (with ni^cr ~ cl ^Ci^^): 

which constrain the hopping of electrons (with hopping amplitude Uj) to the low energy sector of the Fock space, 
i.e. doubly occupied sites are not created. In the language of radial-gauge slave bosons we have 

St^J2^^ i»j n,n+l/./,„+i,^/j,n,o-rj,„+l (26) 
n i.j,(T 

using the relations Eq. (p^). A formal solution of a lattice model with arbitrary hopping is possible since the action 

is bilinear in fa'' and it is found in analogy to the atomic limit of the previous section. The ^—derivatives, which 
enforce the constraint, are now taken at each time step and each lattice site: 

(27) 

with the fermionic inverse propagator matrix 

Dij^nn'[£,i,n\ = ^iJ-^Lrin' + ^ tijf j,n+iri^n+lSn,n' + l (28) 
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and its local part: 



DfX.' = Kw - ^^,ne^''L,n' + l (29) 

It is understood that the anti-periodicity is taken care of implicitly. Although the sum over — 0, runs over 
all possible values of "hole occupancies" on each site and at each time step, "unphysical paths" will be projected out 
due to the product of ^—derivatives. Hereby, unphysical paths are the disconnected paths of the variable r,; „ 
This observation simplifies the calculation considerably. An explicit example is presented in Appendix where the 
partition sum of the two-site cluster is calculated. 

Secondly we consider non-local interactions. The infinite-C/ Hubbard model, that we considered so far, is a highly 
idealized model which is unlikely to be realized in nature. In the attempt to set up a lattice model for interacting 
electrons, one has to realize that many other non-local interaction terms may contribute, (for a recent discussion see 
VoUhardt et al. p^]). On top of the expected density-density interactions, we also encounter spin-exchange terms, 
correlated hopping terms, four-site terms, etc.. Here, using our representation, it turns out that the density-density 
interaction term, 5nioc — X^n Si j a a' ^ '^i,n,cr n-j,n,tT' J can be included without introducing a fermionic interaction 
term. Rather, the contribution to tke action may be rewritten in the form: 

n i,j 

where we used the identity = ri^„, valid in the physical subspace. This Gaussian form of the "interaction term" 
has to be contrasted with the canonical fermionic representation and the Cartesian representation for slave bosons 
where the action remains a true interaction term. In both latter cases, it cannot be treated in a straightforward 
fashion. However, here it appears as additional Boltzmann factors in the expansion of the partition sum Eq. 

The last detail that prevents us from taking advantage of both the radial representation and a conventiona 
agator expansion is the integration range for the bosonic field. Indeed it would be very tempting to extend it from 
[0, oo] to [— oo, oo], such that one could work in momentum space. However such a step is unlikely to yield a mean- 



prop- 



ingful answer, since double occupancy is related to the value = — 1 (see Eq. (15)), which would then be inside 
the integration range. This difficulty may nevertheless be circumvented by i) introducing a new integration variable 
Xn = r'^, ii) making use of the fact that — r„ in the physical subspace, iii) adding a potential term Wxn{xn — 1) 
that precisely vanishes in the physical subspace, iv) extending the integration range for a;„ to [— oo, oo] and, v) sending 
W to infinity at the end of the calculation that would definitely project out any occupancy different from zero and 
one. This transforms the bosonic part of the action Sh in Eq. (pTj) into: 



Si,=6Y,{{^^n + \o){Xr,~l) + WXn{Xn-l)) ■ (31) 
n 

We skipped e here and in the following equations since x„ now extends to negative values. Integration over Xn with 
measure one, now from — oo to -|-oo, yields, in combination with Eq. (p^, the partition sum in the atomic limit: 



( fr r ^ &(.A„+Ao)^.V.(-±^)^^ 



N 



-I 2 



1+n 



g-(5(iA^+Ao-/A) 



rn — 1 



= hm fl + 2e^^ + e2'3(A'->^)') (32) 

As anticipated above double-occupancy is annihilated by taking the limit — > oo while, in the opposite limit — > 0, 
i.e. hadn't we included this potential term as is customary in slave boson calculations [ p3[ , the third contribution 
could well carry a substantial weight |34|. This scheme can be extended to the calculation of the Green's function, 
and in contrast to the above calculation (Eq. (|3^)), W does not enter the result for Zi,tGat,a{T), as explained in 
Appendix Consequently, in a hopping expansion of the partition sum, W would be related only to those paths 
where there is no change in occupancy on one or several sites over the entire imaginary time interval. 

Finally we apply this scheme to lattice models. The measure of the functional integral is the same as for the single 
site problem except that all fields now carry a site index. The potential term, introduced above to enforce physical 
values for the amplitude of the slave bosons, may be implemented as a global term. It is independent of the site index. 
In summary, we find for spin-i fermions on a lattice, interacting through an arbitrary Coulomb-like interaction, which 
is locally infinitely strong, the following expressions for the action in the language of radial-gauge slave bosons: 

S^Sf + Sb + St (33) 
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where: 



i.n,a 

-Sfc = {iK^n{xi,n - 1) + W/'a;i,„(a;j,„ - 1) + - Xi^„)(l - Xj,n)) 

i.n 



(34) 



n i,j,cr 

and the partition sum is given by: 



2tt 



dxi 



(35) 



The measure is now trivial, and the interaction terms included in S^, are bilinear. As an example of how this method 
works for a lattice problem we solve the Ising chain in Appendix 



V. CONCLUSIONS 

Models of strongly interacting lattice fermions have been the focus of many publications on the Mott-Hubbard 
metal-insulator transition, itinerant magnetism and strange metallic phases in recent years. Apart from the thorough 
analysis of the Anderson impurity case |^^, slave boson approaches have been mainly restricted to saddle point 
evaluations. However low energy spin fluctuations, for example, cannot be implemented adequately, as quantum spin 
fluctuations are only included as perturbation. Furthermore, this technique is not suitable to discuss non-Fermi liquid 
behavior, as corrections of higher order in 1/M have to be included. 

This necessitates to further pursue non-perturbative slave boson approaches. Here we considered the construction 
of a functional-integral formalism in the radial gauge which is not restricted to fluctuations around the M — oo 
saddle point. It is defined on a discretized time mesh, and has a well-defined continuum limit. We restricted our 
considerations to models in the limit of large (positive) on-site energy, that is, we work in the constrained Fock space 
with no double occupancy. 

These models allow, in the radial gauge, to perform the functional integrals with real bosonic fields. The phase(s) 
can be absorbed in the constraint Lagrange multipliers which are thereby promoted to real (time-dependent) fields. 
Dynamics of the radial fields themselves does not exist (as discussed in Appendix ^). They serve to enforce the 
constraint for each time step by keeping track of the motion of empty sites, as exemplified for the local problem in 
Eqs. ( p^ ) and (p3|), and for the itinerant problem in Eq. (p8|). This procedure can be extended to other slave boson 
representations as, for example, to the Kotliar-Ruckenstein slave bosons, the representation of which is given for the 
Ising chain in Appendix ^ . 

Non-local terms in general render the functional integral unsolvable. Yet Coulomb-like terms with non-local density- 
density interactions can be rewritten in bilinear form with radial fields, that is, they represent Gaussian terms in the 
functional integral (Eq. ( |30| ) in Section IV, and also other non-local interactions may be rewritten in bilinear form, 
e.g., the Ising spin coupling within the Kotliar-Ruckenstein representation, see Appendix This is neither possible 
in the canonical fermionic representation nor in the Cartesian representation for slave bosons. These bilinear terms in 
the Lagrangian generate a finite dispersion for the (radial) slave particles. In order to keep the bosonic part Gaussian 
and to allow for a momentum space representation, one has to extend the range of integration for the square of the 
radial field to minus infinity. It can be implemented with the observation that the functional integration stays correct 
if we add an additional global Gaussian term, that is, a potential term W{xi^n — l)a^i,Ti where Xi^n is the square of the 
radial part of the slave boson, and send ^ oo at the end. This proved to be a valid procedure to constrain the 
amplitude of the slave boson to the physical values and 1, without introducing any additional complications (see 
Appendices |^ and ^ and the paragraph with Eqs. ( ^ and (^) in Section tV). 

In this article we advanced two mechanisms for the solution of a functional integral for constrained electronic 
problems: either we keep continuous fields which are originally introduced on a discrete time mesh but have a well- 
defined continuum limit, or we enforce the constraint in the first step and thereby reduce the slave boson radial field to 
an Ising- like variable in a discrete representation. The first approach is summarized in the final expressions Eqs. (jsj) 
and (^5|) — and for Kotliar-Ruckenstein slave bosons in Eq. (C2) — in which the bosonic fields are real fields and 
the bosonic action is of bilinear form. The second approach is discussed in Section III, and extended to itinerant 
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problems in Section IV, Eqs. ( p7[ ) and (|29|), which can be taken as a starting point for a numerical evaluation in the 
strong coupling limit that does not build on a Hirsch decoupling. Furthermore, non-local interactions can be included 
without the need of any additional decoupling. 
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APPENDIX A: 

Below we motivate why the radial amplitude field has no dynamics and, furthermore, why the limit of integration 
is shifted to the negative real axis by an infinitesimal amount e. We start from Sb Eq. in Cartesian representation. 

n 

and substitute &„ — r„e''''". 

5*6 - ^ [r„(r„ - r„_i) + r„r„_i (l - e'*"""'* (1 - iSX)^ 

n 

= i6{ipn + X)[rl ~ Srnf-n] + Sr„r„ + "further terms in 0(6^)" + 0{S^) (A2) 

n 

In the last line we observe that J2n ^'''n'f'n — > f '^'^^ ~ ^ ("^^^ periodic boundary conditions for the bosonic 
field). The "further terms in 0(<5^)" can all be neglected with respect to the linear terms. Here we introduced the 
time-dependent constraint field A„ = -I- A. The included term quadratic in b is special: if ~8rnTn is positive, it 
is always a negligible correction to r^^. However, for Srnrn > it may interpreted as an infinitesimal shift e of to 
negative values for r„ — > 0. Consequently, the leading term in St now reads: 

5b = ^i<5A„(r2 -e) . (A3) 



There is no dynamic term of the field r„ — as there is one for &„ in Eq. ( Al ) — since this amplitude field is a real field. 
On account of the constraint we avoid the essential difficulty that one encounters when one uses polar coordinates as 
discussed by Edwards and Gulyaev ||30[| . Particularly, the part of the action in polar coordinates, that corresponds to 
the Gaussian part of the action in Cartesian coordinates, cannot be written in a Gaussian form. This is true even to 
first order in S, and in Jsot higher order terms in the polar fields must be kept as well. 

APPENDIX B: 

As an illustration of how our discrete representation yields the partition sum we explicitly calculate it for the two- 
site cluster. An instructive "warm up" exercise consists in discretizing the imaginary time axis into six time steps. In 
order to evaluate Z Eq. (|27| ) we first choose the physical path: 

ri = r2 = 1; ^3 = r4 = 0; = rg = 1 on site 1 

fi = 0; f2 = . . . = r5 = 1; = on site 2 (Bl) 



10 



Here and in the following we use i„ = Xi=2,n for anything on site 2. To built up the inverse propagator matrix we 
observe that all ^'s for which r is equal to 1 can be set to zero, since the derivatives with respect to those ^'s have to 
be taken from the prefactor. As a result the inverse propagator matrix is: 



D = 



1 
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-5t 






















1 





































1 






































1 














-5t 























1 






































1 
















-5t 




















1 






































1 






































1 

















-5t 




















1 





































1 / 



(B2) 



Along this path the hole hops from site 1 to site 2 at time 2, and back to site 1 at time 5. Accordingly one would 
expect that t only enters D on two entries: (2, 1) and (4, 5). But it also appears at two other entries. Nevertheless the 
latter two entries, which do not correspond to physical processes, are easily seen not to contribute to the determinant. 
The latter reads: 



det [D] = l + iStfe^'^^3^i^i^e 



(B3) 



If we now let the number of time steps to be TV and gather all second order processes in t (their number being 
N{N — 1)), we obtain their contribution to the partition sum as 



Z2 = 2 {(3tf e(~-2)^^ 



(B4) 



as it should. With some patience one may extend this procedure to an arbitrary number of hopping processes, and 
we checked that they can all be summed up to 



Z = l + 2e''^2cosh(/3i) 



(B5) 



which is the correct answer. 



APPENDIX C: 

In this appendix we illustrate the method described in Section for the Ising chain. To that aim we make use of 
the Kotliar and Ruckenstein representation of the Hubbard model Q in the limit L/ — > cxd. We thus introduce two 
auxiliary fermions /| and /|, and three auxiliary bosons e, representing empty sites, and and pi representing singly 
occupied sites. On each site i they are subject to three constraints: 

p{aP^■.^ = /L/- <y -t, i (CI) 

which are respectively enforced by three Lagrange multipliers denoted by ai and Xi^a. The phase of all three bosons 



can be gauged away [^,^ 37|, and the three Lagrange multipliers are promoted to fields. Since, for the Ising model, 
we are working at exactly one electron per site, the e-field, representing empty sites, has to be fixed to zero. We also 
could implement this "half-filling condition" in the standard way by introducing a chemical potential for the electrons. 
However here, it is more convenient to introduce a Kronecker (5- function q in the measure (see Eq. (|C3| ) below). 
Just setting the e-field to zero from the outset would introduce a spurious divergence in the measure. The action is 
thus: 
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Pi,n.(T ^i.7i 1 I 

i,n,(T i,n \ (7 / 



(C2) 



It consists of a fermionic part, the first contribution, and the remaining ones form the bosonic part Sb- Here, Jij is 
the Ising-spin couphng, and W was introduced in Eq. (^) and the preceding discussion. After the integration over 
the fermions we obtain the partition sum as: 



z = hm <^ n 



N 

n 

,n=l 
N 



Sdai 



2tt 



n 



5d\i 



2tt 



ei,„,0 



n 1+n 



-i5\i 



(C3) 



We first integrate over the constraint fields Ao-, then over a and finally over e. As a result the term with W drops out — 
since the Ao--integrations already enforce the constraint to the point that Pi^n,a- can take only the discrete values and 1. 
Furthermore we observe that pi_n,a loses its time-dependence since the previous integrations resulted in the two straight 

"world lines" n„ HPin^)S{pini ^ 1) + Tin HPin^ - l)<5(Pmi) for each site i. Thus, setting pi^^ = p^^i^^ = ... = Pi^^a 
we get: 



n 



dp. 



dPi.i {^&{PtA)HPtA - 1) + HPtA - 



,- Ji,i(.Pi.-[-Pi,l)iPj.-^-Pj.l) 



(C4) 



which holds for any coupling matrix Ji j and topology. If we now restrict ourselves to nearest neighbor interaction 
and a chain of length L with open boundary conditions, we can recursively integrate over Pi,x, • • ■ , PL.^^ and 
PL, I to obtain the known result 



Z = 2 (2 cosh {(3 J)) 



(C5) 



APPENDIX D: 



We calculate the Green's function in the atomic limit using the action Eq. (^4|) in which a global constraint was 
introduced so that the radial amplitudes can be integrated from — oo to oo. In combination with the expression of 
the Green's function Eq. (|20|) and the propagator Eq. (p^) we get: 



Za.tGa.t.aim 



[No 

1) = hm n / 



dXr 



5dX 



27r 



!ig-5(iA„(.„-l)+T^.„(.„-l)) ^^^^^^ 



(Dl) 



This relation is easily verified either by straightforward integration over the Grassmann fields or, more directly, by 



realizing that in Eq. (|2l|) one has to replace ^„ 



i5X„ 



in the propagator D and in detl?, and is 



replaced by the bosonic path integral e '^('-^'•(^" i)+Wx„{x„ i))^ Integration over A yields: 

^atGat,.(™-l)= lim e*>('"-iMrr / dx„e-^^^"(^"-i) 

X ^(xi - 1) n ^(^") n ^(^« - 1) + e^^^(xi) n '^(^" + 1) n -^(^ 



\n—'m-\-l 



\n—m-\-l 



— g'5/x(m-l) 



(D2) 
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The second contribution above, which is related to double occupancy, cancels due to the additional projection in 
the external electron operators, represented by the factors xiXm+i in Eq. (^o|). For this reason W drops out of the 
calculation of the Green's function, in contrast to the calculation of the partition sum. This holds for the calculation 
of higher correlation functions as well. 
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Figure captions 

Figure 1: Contributions to tlic atomic partition sum for infinite on-site interaction U . 

Figure 2: Time evolution of the fields r and /o- in the imaginary-time local Green's function, Ga{m — 1). The 
dashed line is a "particle chain" (see text for details). "Hole chains" are not depicted. 



14 



" Permanent address: Laboratoire Crismat-ISMRA, 6, Bid. du Marechal Juin, 14050 Caen Cedex 4, France. 
[1] Y. Vilk, and A.-M. S. Tremblay, J. Phys. I France 7, 1309 (1997). 

[2] P. W. Anderson, The Theory of Superconductivity in the High-Tc Cuprates (Princeton University Press, Princeton, 1997) . 
[3] S. Pairault, D. Senechal, and A.-M. S. Tremblay, Phys. Rev. Lett. 80, 5389 (1998), and condmat/9905242. 
[4] Ph. Brune and A. P. Kampf, condmat/0001210. 

[5] In a different approach to these effective models, a degenerate perturbation theory is set up ab initio so that a Hamilton 
operator in the projected Fock space, restricted to low-energy excitations, can be determined: M. Takahashi, J. Phys. C 
10, 1289 (1977). 

[6] E. O. Tiingler, and T. Kopp, Nucl. Phys. B 443, 516 (1995). 

[7] S. E. Barnes, J. Phys. F 6, 1375 (1976); F 7, 2631 (1977). 

[8] G. Kotliar, and A. E. Ruckenstein, Phys. Rev. Lett. 57, 1362 (1986). 

[9] R. Fresard, and P. Wolfle, Proceedings of the Adriatico Research Conference and Miniworkshop "Strongly Correlated 

Electron Systems III", Eds. G. Baskaran, A. E. Ruckenstein, E. Tosatti et Yu Lu, Int. J. of Mod. Phys. B 6, 685 (1992). 
[10] However Bose condensation may be suppressed, even on a mean field level, and bond pairing of bosons emerges instead if 

the constraint is recast into an interaction term between the various slave particles: T. Kopp, F. J. Seco, S. Schiller, and 

P. Wolfle, Phys. Rev. B 38, 11835 (1988). 
[11] W. Metzner, and D. VoUhardt, Phys. Rev. Lett. 62, 324 (1989); W. Metzner, Z. Phys. B 77 253 (1989); W. Metzner, and 

D. VoUhardt, Phys. Rev. B 37, 7382 (1988). 
[12] L. Lilly, A. Muramatsu, and W. Hanke, Phys. Rev. Lett. 65, 1379 (1990). 
[13] R. Fresard, M. Dzierzawa, and P. Wolfle, Europhys. Lett. 15, 325 (1991). 
[14] W. Zimmermann, R. Fresard, and P. Wolfle, Phys. Rev. B 56, 10097 (1997). 
[15] A. Tandon, Z. Wang, and G. KotUar, Phys. Rev. Lett. 83, 2046 (1999). 
[16] G. Seibold, E. Sigmund, and V. Hizhnyakov, Phys. Rev. B 57, 6937 (1998). 
[17] R. Fresard and G. Kotliar, Phys. Rev. B 56, 1290 9 (1997) 



[18] H. Hasegawa, J. Phys. Soc. Jpn. 66, 1391 (1997); |gond-mat/0005271 . 
[19] M. Lavagna, Phys. Rev. B 41, 142 (1990). 

[20] Y. Okimoto, T. Katsufuji, Y. Okada, T. Arima, and Y. Tokura, Phys. Rev. B 51, 9581 (1995). 

[21] P. Coleman, Phys. Rev. B 35, 5072 (1987). 

[22] N. Read and D. M. Newns, J. Phys. C 16, 3273 (1983). 

[23] E. Arrigoni, C. Castellani, M. GriUi, R. Raimondi, and G. C. Strinati, Phys. Rep. 241, 291 (1994). 
[24] L. Gehlhoff, and R. Fresard, J. Phys: Condens. Matter, 8, L13 (1996). 

[25] J. Kroha, P. J. Hirschfeld, K. A. Muttalib, and P. Wolfle, Solid State Commun. 83, 1003, (1992); J. Kroha, P. Wolfle, and 

T. A. Costi, Phys. Rev. Lett. 79, 261 (1997); J. Kroha, and P. Wolfle, Acta Phys. Pol. B 29 (12), 3781 (1998). 
[26] J. E. Hirsch, and R. M. Fye, Phys. Rev. Lett. 56, 2521 (1986). 
[27] N. E. Bickers, Rev. Mod. Phys. 59 845 (1987). 
[28] N. E. Bickers, Ph.D. thesis, (Cornefl University) (1986). 
[29] P. W. Anderson, Phys. Rev. 110 827(1958); Phys. Rev. 112 1900(1958). 

[30] S. F. Edwards, and Y. V. Gulyaev, Proc. Royal Phys. Soc. A 279, 229 (1964); K. Ito, Mem. Amer. Math. Soc 4, 1951. 
[31] P. Coleman, Phys. Rev. B 29, 3035 (1984). 

[32] A connected path is a temporal succession of r-values in which an initial value of r equal to 1 at a site is continuously 
"transported" through the lattice until the final time step. For example, a "transport process" for which the value of rj at 
site j jumps from at time step n to 1 at time step (n -|- 1) is valid only when it is connected to — via a matrix element ti,j 
— a site i where the value of jumps from 1 at time step (n -I- 1) to at time step {n + 2). For disconnected paths such 
jumps are not related, and they are excluded since such paths cannot generate "complete chains" (see previous section). 

[33] D. VoUhardt, N. Bliimer, K. Held, M. KoUar, J. Schlipf, M. Ulmke, and J. Wahle, Z. Phys. B 103, 283 (1997). 

[34] For the atomic limit, this third term in Eq. (|^ with ^ in fact is equal to the contribution from the state with double 
occupancy for zero on-site interaction energy. However, this observation does not imply that we would have the option to 
recover the free limit without on-site interactions for an itinerant model by applying W ~> 0. The phase relation between 
the various local states would still be missing. Therefore, the third term in Eq. (p2|) is spurious and we always have to take 
the valid limit — > oo. 

[35] The original Kotliar and Ruckenstein representation has been introduced for problems with finite on-site interaction. In 
that case an additional field representing double occupancy needs to be introduced. Direct, but lengthy, calculations show 
that this field needs to be complex in order to get the correct result for the partition sum and correlation functions in the 
atomic limit. 

[36] Th. Johcoeur, and J. C. Le GuUlou, Phys. Rev. B 44, 2403 (1991). 



15 



[37] Y. Bang, C. CastcUani, M. Grilli, G. Kotliar, R. Raimondi, and Z. Wang, Proceedings of the Adriatico Research Conference 
and Miniworkshop "Strongly Correlated Electron Systems III", Eds. G. Baskaran, A. E. Ruckenstein, E. Tosatti et Yu Lu, 
Int. J. of Mod. Phys. B 6, 531 (1992). 



16 



Fig. 1 




Fig. 2 




